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Abstract. Phonon lineshapes in atom-surface scattering are obtained from a simple 
stochastic model based on the so-called Caldeira-Leggett Hamiltonian. In this single- 
bath model, the excited phonon resulting from a creation or annihilation event is 
coupled to a thermal bath consisting of an infinite number of harmonic oscillators, 
namely the bath phonons. The diagonalization of the corresponding Hamiltonian 
leads to a renormalization of the phonon frequencies in terms of the phonon friction or 
damping coefficient. Moreover, when there are adsorbates on the surface, this single- 
bath model can be extended to a two-bath model accounting for the effect induced by 
the adsorbates on the phonon lineshapes as well as their corresponding lineshapes. 
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1. Introduction 

Bulk and surface phonon dynamics is a very active field of research. Usually the 
information about this dynamics arises from scattering experiments, where the typical 
observable quantity is the lineshape associated with a surface excitation, i.e. a creation 
or an annihilation event. At present, the available experimental techniques, such as the 
surface scattering with He-atom beams [1] or the He spin echo [2], allow to explore higher 
and higher angular and energy resolutions. Therefore, it is important and necessary 
to have at our disposal appropriate theoretical formalisms which allow us to extract 
reliable information from the experiment, such as phonon lifetimes and frequency shifts, 
phonon-dispersion relations or multiphonon backgrounds, for example. In this regard, 
a qualitative explanation was formerly provided for multiphonon scattering processes 
in terms of a semiclassical formula for energy and momentum exchange within the so- 
called trajectory approximation as a generalization of the Brako-Newns formula [3,4]. 
Levi and Bortolani [5] also developed a general multiphonon formalism based on the 
time-evolution approach. More recently, Gumhalter [6] has made a revision on single 
and multiphonon atom-surface scattering within the quantum regime. However, as far 
as we know, little attention has been paid to the theory of phonon lineshapes. 

In this work we develop a theory of phonon lineshapes within a stochastic approach 
based on the so-called Caldeira-Legget Hamiltonian [7], which constitutes the paradigm 
of stochastic dynamics and system-plus-reservoir approaches. In particular, we focus 
on an excited phonon in the presence of a phonon field and study its relaxation. To 
this end, in next Section we shortly review the theory of Manson and Celli, based on 
the transition matrix formalism [8-10], in order to introduce the notation and their 
main theoretical findings on single- and multiphonon scattering. In Section [3J the 
Caldeira-Leggett Hamiltonian is introduced to develop a theory for damping phonons 
within a single-bath model (phonon bath). This Hamiltonian, which displays a high 
degree of versatility, has been used by our group to describe surface diffusion of single 
adsorbates [11-13] and interacting adsorbates within a two-bath model (phonon and 
adsorbate baths) [14,15]. This Hamiltonian has been applied, also recently, to classical 
atom-surface scattering to describe angular distributions [16,17] and energy loss spectra 
[18,19]. Finally, in Section @J different aspects of this new approach are discussed and 
applied to He scattering off a Cu(001) surface with a 3% CO coverage at low (surface) 
temperatures [20]. 

2. A brief account on the transition matrix formalism 

The observable in atom-surface inelastic scattering experiments is usually the so-called 
differential reflection coefficient, which gives the fraction of probe particles (e.g. He 
atoms) scattered off into a final solid angle dQf and an energy interval alEf. At a 
theoretical level, this observable is obtained as a transition rate, dividing the incident 
flux crossing a plane parallel to the surface and then multiplying by the density of 
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available final states of the probe particle [8-10]. This allows us to express the transition 
rate as [8,9] 

i r+oo 

Tl(k f , ki) = ± / e -^-^ )t/h \r{k f , k t )\ 2 ^( e -k.[R i+ u ;( o)] e ,k.[R J+ u jW]) rf£) (1) 

h J ~°° 3,1 

where e« and ej are, respectively, the initial and final particle energies, tko = ej — e« is 
the energy exchanged in the scattering process, k = kf — kj is the total wave vector 
exchanged by the scattering particles, and rjj = r(k/,kj) is a pairwise transition 
operator. The (time-independent) equilibrium position of the jth surface unit cell 
is given by Rj, with j being a discrete two-dimensional variable for each layer. For 
simplicity, we will assume that there is only one atom per unit cell and, therefore, Uj 
corresponds to the displacement of the jth atom with respect to Rj due to its vibrational 
motion. In (pQ), the sums extend over all Rj and the average is over all initial lattice 
states (once the sum over all possible final lattice states has been carried out). 

Equation ([TJ) has the general form of the fundamental equation for a large class of 
scattering treatments involving a many-body target, 

n(k f ,k { ) = ^\r(k f ,k l )\ 2 S(k,u). (2) 

This equation is the product of a form factor, |r(k/,k,)| 2 , for the scattering centers 
and a dynamical structure factor, S(k, to), which depends on the average over scattering 
center positions. The latter factor describes, precisely, the lineshapes in the kind of 
experiments we are interested in here, and can be expressed as 



S(k,u) = J I(k,t)e wt dt, (3) 

i.e. as a time Fourier transform of the so-called intermediate scattering function, 

I(k, t) = J2 e- ik - Ri e ik - R ^' ( e -*-ui(o) e ik. Uj (t)^ ( 4 ) 

3,1 

If the harmonic approximation is valid for any vibrational lattice mode, then the 
(vibrational) average in (Hj) can be easily carried out by standard methods, which yields 

( e -ik-Ui(0) e ik-Uj(t)^ = e -T !f (k,^ ( 5 ) 

where 

Ty(k, t) = 2W(k) - ([k • u,(0)][k • Uj(t)]>, (6) 

and 

iy(k) = i([k- Uj (t)] 2 ) = i([k- u ,(0)] 2 ). (7) 

The factor W(k) is the argument of the so-called Debye-Waller (DW) factor, e~ 2W ^ k \ 
which governs the decay of the scattered intensity as a function of the surface 
temperature through the canonical ensemble average. This contribution to the thermal 
attenuation of the scattered intensities arises mainly from the momentum transfer 
perpendicular to the surface. 
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The second term in ([H]) accounts for the energy transfer between the surface and 
the scattered particles, which is governed by the displacement autocorrelation function 

^(k,t) = ([k.u,(0)][k. Ui (t)]}. (8) 

Here we are dealing with the problem of inelastic scattering and, therefore, our interest 
relies on the evaluation of this displacement correlation function. Hence we express the 
vibrational amplitude of the jth scattering center in terms of surface normal modes as 

= £ A Qi „e(Q, u) [a Q At)e- lQ ^ + al^)^ ^] , (9) 



where 



and 



OQ,u(t) = a Q , u (0)e~ lu;AQ)t , (10) 



In these expressions, Q is the parallel-phonon momentum, v represents additional 
phonon quantum numbers (e.g. the Rayleigh-mode branch or normal momentum of 
bulk modes), M is the mass of the surface atom in the unit cell, e(Q, v) is the 
polarization vector, and a^ )U {t) and <Y_Q U (t) are the time-dependent phonon annihilation 
and creation operators, respectively. The polarization vector is related to the frequency 
distribution function, p a/3 , as [8,9] 

Pa/3 = £ e Q (Q, v) ■ e /3 (Q, u)S(u- w„(Q)), (12) 

where a and /3 are the Cartesian components of the polarization vectors. From ([9]) and 
f fTOj) . and after performing the corresponding thermal averages, OH]) becomes 

lfy(M) = ^A^[k.e(Q,z/)] 2 {n l/ (Q)e^( R '- R J )e^^ t 

+ K(Q) + i] e -iQ-(^-R-i) e ^(Q)*| ; (13) 

where ^(Q) = [ e ^(Q)/ fc s T _ ij- 1 i s the Bose-Einstein factor. By Taylor expanding the 
exponential of £7y(k, t), with Ty decomposed as in the multiphonon contributions 
associated with the total inelastic scattering process can be obtained. The first term 
of this series expansion (namely the unity) gives the diffraction intensities, i.e. the 
elastic contribution to the scattering process, after carrying out the integral (j3]) and 
then computing the transition rate (|2]). At this level of approximation, as reported by 
Manson [9], these diffraction intensities will be given by the product of three factors: the 
DW factor, the form factor and the structure factor. In particular, the (static) structure 
factor is the time-independent factor arising from (pEJ) at zero temperature, 

5(k) = £ e - ik - R ^e ik - R! , (14) 

U 
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with the diffraction peaks then expressed in term of 5-functions. The next term, 
the linear one, provides the single-phonon contribution, with its corresponding rate 
containing an extra factor, namely the Bose-Einstein factor. Also in this case inelastic 
intensities are given by 5-peaks. Finally, the multiphonon contribution is obtained from 
the remaining terms of the Taylor expansion or, equivalently, by subtraction of the zero 
and single-phonon contributions. 

Since we are interested in phonon lineshapes, a damping mechanism has to be 
introduced in a theoretical model in order to replace the 5-peaks by Lorentzian 
lineshapes. The goal in next Section is to assume a simple model, the so-called Caldeira- 
Leggett (single-bath) model. The corresponding Hamiltonian [7] has been widely used 
and applied in the literature to a large variety of physical processes in presence of 
dissipation [21]. More specifically, this approach has been extensively used in the 
so-called Kramers' turnover problem [22, 23] describing the escape from a metastable 
potential, in the vibrational dephasing problem of small molecules in liquids [24], and 
in diffusion problems [11-15] with no application to lineshapes. 

3. A simple model for phonon lineshapes: The Caldeira-Leggett model 

Consider the damping of a single phonon (system) excited by some external source 
and coupled to a thermal bath at a given temperature (heat bath). This dissipation 
mechanism is not very sensitive to the particular loss process involved and, therefore, the 
heat bath can be assumed as consisting of an infinite number of independent oscillators 
(phonon field) in thermal equilibrium at the surface temperature. Thermal fluctuations 
in the bath or reservoir, which feed with noise the excited phonon, are related to the 
friction coefficient or damping constant through the well-known fluctuation-dissipation 
theorem. The system is usually represented by a harmonic oscillator weakly coupled 
to the reservoir and undergoing an energy flow towards the reservoir. Conversely, 
the reservoir fluctuations also couple back into the system. The driving force (He 
atoms) exciting the phonon is assumed to be a linear perturbation and, at zero 
order, no influence in the subsequent phonon dynamics is expected. From the linear 
response theory, the dynamical susceptibility will be that of a particle subject to a 
one-dimensional harmonic potential [14,15]. 

Taking into account the comments above, the full Hamiltonian describing the 
process can be written as [7,21] 



where Hs describes the Hamiltonian accounting for the system dynamics (i.e. the excited 
phonon), Hr is the Hamiltonian associated with the thermal bath and H$r represents 
the coupling between system and bath. The so-called Caldeira-Leggett Hamiltonian has 
the structure given by (fl5l) and is usually expressed in the form 



H = H s + Hr + H S r, 



(15) 




(16) 
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where (q, p) and m represent the system variables (coordinate and conjugate momentum) 
and its mass, respectively; V(q) is the external potential acting on the uncoupled system, 
which here we assume to be harmonic; (xi,pi) stand for the variables of the zth bath 
oscillator, characterized by a mass m ; and a frequency uf, and q is the coupling constant 
between the ith bath oscillator and the system. 

Within this single-bath model, tracing over the bath degrees of freedom in the 
equations of motion (expressed in the Heisenberg picture) makes evident the damped 
motion undergone by the system coordinate [21]. This tracing gives rise to a quantum 
generalized Langevin equation, 

mq + ^l + m [ 7 (t - f) q{t') dt' = £(t), (17) 



dq 

which describes the evolution of the system coordinate (q) under the influence of the 
external harmonic potential V(q), a friction kernel j(t) and a random force £(i) coming 
from the bath thermal motion. In this model, the friction kernel (phonon friction) is 



1 N c 2 
nf(t) = — ) — % —~ cos Uit, 

i=i 1 



while the spectral density characterizing the bath oscillators is 

N 2 

"M=*£^%- <)■ (19) 

i=l 

On the other hand, the random force reads as 

m = E c i { - ^ ?(«)] cos ^t + 1® sin ^4 > (2°) 

j=\ ^3 J 3 3 J 

which satisfies the conditions for a Gaussian white noise provided the friction is Ohmic, 
i.e. 7(£) = 27<5(t). In this case, Eq. (TlTl) becomes a standard quantum Langevin equation, 

dV(q) 

mq + -— ^ + m 7 g(t) = f(t). (21) 
oq 

Quantum-mechanically, the Hamiltonian ffl6|) for a harmonic potential V(q) = 
mu 2 q 2 /2 can be expressed in terms of system (a, a^) and bath (bj,bj) operators as 

N N 

H = frw G (a)a + 1/2) + ^ ^j{b)bj + 1/2) + h^K^ab] + atbj + a)b] + abj) 

j=i j=i 

n N c 2 

+ -^—{a^ 2 + a 2 + 2a^a + 1) V — ^, (22) 



with 



c 



kj = — 3 (23) 
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accounting for the strength of the bilinear coupling between the system and the bath 
oscillators. As is apparent, two quantum transitions in the system are allowed according 
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to the existing second order terms. The Hamiltonian ( 12 2 p can be diagonalized via 
a normal mode transformation with mass weighted coordinates. The corresponding 
normal- mode Hamiltonian can be then written as a sum of iV+1 oscillators [22-24], 

B = \P 2 + \ <Q 2 + t( 1 2 P !+ 1 2 »M) • P") 

3=1 V J 

where (w ,Wj) denote the new system and bath frequencies (eigenvalues) and (Q,P) 
and (Qj,Pj) are the new system and bath variables, respectively. For Ohmic friction, 
the explicit expression of w is 



w = ±y / u; 2 - 7 2 /4 + i 7o /2 = ±w + *7o/2. (25) 
Similarly, for Wj, 

Wj = ± J to] - 7 |/4 + z 7i /2 = ±Wj + i 7j -/2 . (26) 

The real parts in both cases give renormalized frequencies, Wq and Wj, and the 
imaginary ones the damping constant or phonon friction coefficient. Notice that the 
new frequencies are no longer a simple linear function of the friction coefficient, as it 
usually happens when only single quantum transitions are allowed. Furthermore, after 
diagonalization, we have a set of iV+1 damping phonons. This procedure is equivalent 
to write Hamiltonian ( 1241 in terms of real frequencies Wk [25,26]. Hence, in terms of 
operators, the Hamiltonian (|24|) can be recast as 

N 

H = J2^M4a k + l/2), (27) 

k=0 

where (a^,a) are related to the new phonon variables (Qk,Pk)- The time evolution of 
such operators is given by 

a k {t) = a k (0)e iw " t , (28) 

with the minus sign usually chosen for the real part of the frequencies in ( 1251) and ( 1261 . 

In order to obtain the phonon lineshapes, we proceed as previously shown, although 
now the vibrational amplitude of the jth scattering center is given in terms of the new 
normal modes or damping phonons, expressed through (a\,ak) as 

= E AQ,,e(Q, u) e-^' 2 [a Q At)e- iQ ^ + al Q<u (t)e^] , (29) 

where, for simplicity, we keep the same notation (Q, v) for quantum numbers, as in (jHJ), 
and 



Aq '-^2(JV + 1K(Q)- (30) 

Note that now the time-dependent operators are thus given by (128]) and the 
corresponding frequencies by (1251) and (|26|) . unlike Eq. (TlOl) . where frequencies were 
real. In this way, acquires formally the same expression, 

Ty(k,t) = ^([k • u,(0)] 2 ) + i([k • u,(t)] 2 ) - MM), (31) 
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but with the displacements being given by (1291) . This fact has important consequences: 
though the first term in (13T1) is W, as in ([7]), the second one is no longer independent 
of time (it depends on time as e -7 ^^*) and, therefore, will not contribute as a new W 
term. Regarding the third term, it is also time-dependent, containing the new extra 
factors e -7 "^*/ 2 , e tWu( ^ >t and the complex conjugate of the latter. In particular, this 
displacement autocorrelation function, given above by (fl3"j) . can now be expressed taking 
into account the new normal modes as 



fy(M) = ^A^[k-e(Q,z,)] 2 {^(Q)e iQ ( R - R ^e-^( Q V 



+ K(Q) + i] e -iQ<^i-^) e i^(Q)t e -iAQ)t/2\ (3 2 ) 



i • 

where w u (Q) can be identified by Wk- Observe that this equation is similar to ffTB"]) except 
for the presence of damping factors. Therefore, phonon lineshapes will be obtained from 
the time Fourier transform, as given by ([3]), but considering the damping normal modes 
in terms of (Q, v) and w u (Q,). The interest here relies on the fact that lineshapes can 
be obtained as a function of the new system and bath frequencies issued from a given 
diagonalization due to the presence of the heat bath. This approach thus allows us to 
consider the excited phonon and the phonon bath on equal footing. 

Specifically, the phonon lineshapes obtained from the time Fourier transform ([3]) 
read as 



x , <2(n+m+p) 



v ' ; 2irh A^ll v 2 ' n\m\p\ 1 n 

j,l Q,v n,m,p=0 

x MQ) n MQ) + ip[2n„(Q) + i]V Q -( Ri - R ^ n e - iQ - (Ri - R ^ )m 

(n + m + 2p)7^(Q) 
[u-(m- n)w u (Q)} 2 + [(n + m + 2p) 7l/ (Q)/2] 5 
In principle, this expression contains the main ingredients of the phonon relaxation 
dynamics. For a given phonon, (1331 consists of an infinite sum of weighted Lorentzian 
functions which describes all phonon creation and annihilation events. 



x r. . TZ „n„t-. 77vo2 . i ■ ^37 (r^\ /oi2' 



4. Discussion and applications 



Some consequences of (!33|) can be discussed at different levels: 

(a) Damping phonons in clean surfaces: positions and lifetimes 

The first point to be stressed from (!33j) is the frequency position and lifetime of 
each (creation or annihilation) phonon event. The frequency position is given by f l25l) 
and involves the phonon friction. Thus, the nominal value of the phonon frequency ojq or 
ujj has to be replaced by a renormalized frequency due to the coupling with the phonon 
bath. The phonon lifetime is directly related to such a friction, h/ r y v (Q l ). The different 
excitations of a certain phonon event are given by the sums running over n and m. Each 
elementary contribution is expressed by a weighted Lorentzian shape. Obviously, the 
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most prominent Lorentzian lineshape is related to the most likely phonon event. When 
n = m, we have the diffuse elastic peak. This contribution arises from the creation and 
annihilation of the same number of phonons and is essentially associated with single- 
phonon exchanges (the multiphonon contribution is substantially smaller). The product 
in the quantum numbers (Q, v) shows the contributions of different excited phonons. 
The multiphonon background is given by the different factors appearing in fl33|) . The 
temperature dependence comes through the Bose-Einstein and DW factors. Regarding 
the symmetry of the total lineshape, and according to the detailed balance condition, 
particle energy gains (which requires the system to be in the higher energy state) are less 
frequent than atom energy losses, in agreement with the Boltzmann population factor. 

(b) Damping phonons and the DW factor 

In (133|) . the DW factor appears as e~ w instead of the standard factor e~ 2W . This 
is a consequence of the diagonalization carried out in the Caldeira-Leggett Hamiltonian 
( I16p . The normal modes are now damping phonons given by ( 1281) and, therefore, the 
square thermal averages of the displacements are no longer independent of time. As is 
well known, this factor is a measure of the effect of thermal motions in reducing the 
periodicity of the lattice. It is also a measure of the number of phonons involved in a 
scattering event. In this formalism, it is clearly seen that there is a new contribution 
to the lifetimes (the sum running over p) independent of the creation or annihilation of 
phonons. This contribution comes from the whole damping lattice. 

(c) Adsorbates and phonons 

The diffuse elastic intensity of the He atoms scattered at large angles away from the 
specular direction provides very detailed information on the mobility of adsorbates on 
surfaces. At low coverages, the interaction among adsorbates can be ignored, thus 
allowing to work within the so-called single- adsorbate approximation describing the 
quasielastic and low frequency vibrational lineshapes in terms of the motional narrowing 
effect [27]. The relaxation processes of creation or annihilation adsorbate events are 
due to the coupling of the adsorbate modes (considered as the system) to the phonon 
substrate (single-bath model). Thus, ( 133|) is also valid for describing the relaxation 
dynamics of an excited single adsorbate with a renormalized frequency given by ( T25|) . 
Furthermore, it can be considered as a generalization of previous lineshapes obtained in 
a more phenomenological way [27]. 

Based on the transition matrix formalism, Manson and Celli [8] proposed a quantum 
diffuse inelastic theory for small and intermediate coverages of adsorbates on the surface 
by ignoring multiple scattering effects of He atoms. The dynamical structure factor 
is then obtained by assuming all the lattice vibrational modes (N p h) and point-like 
scattering centers (N ai i) satisfying the harmonic approximation with a given frequency 
distribution function. Therefore, following the same type of arguments, we could assume 
two independent, uncorrected baths to describe diffusion of interacting adsorbates: 
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the first bath consists of the surface phonons and the second bath is formed by N a d 
adsorbates which obviously changes with the surface coverage given by experimental 
conditions [14,15,28]. The same procedure can now be used to obtain lineshapes of 
adsorbates; in particular, the lowest frequency mode or frustrated traslational mode. 
To this end, we have to generalize the Hamiltonian ( Ti~6l) to two baths, one for phonons 
and the other for adsorbates. In this two-bath model (at a given coverage), we take 
one adsorbate as the tagged particle or system, while the remaining ones constitute 
the second bath described by M harmonic oscillators. Thus, when a given adsorbate is 
excited the coupling to phonons and adsorbates simultaneously governs to the relaxation 
dynamics. In this way, the corresponding total Hamiltonian in one dimension will read 
as [14, 15] 



where now the tilde magnitudes refer to the second bath of N a d adsorbates, which are 
also taken as harmonic oscillators. The Cj and coefficients give the coupling strengths 
between the adsorbate (system) and the substrate phonons or other adsorbates, 
respectively. The spectral density for the two baths is defined analogously to the single- 
bath model, 



but now it is split up into two terms, one spectral density associated with the surface 
phonons and another one with the adsorbates. In a similar way, the friction functions are 
defined as in ([TBI , but with the spectral density being (|35|) . The total friction function 
rj{t) also splits into two terms, one due to the phonons, j(t), and another due to the 
collisions with the the adsorbates or collisional friction, A(t): r](t) = ^(t) + X(t) [14,15]. 
After diagonalization of Hamiltonian (13~4"1) and if Ohmic friction is assumed, the new 
renormalized frequencies of the set of N p h + N a d oscillators are expressed in a similar 
way to Eqs. (j25|) and (|26|) but now in terms of r\ instead of 7. A straightforward 
generalization of lineshapes accounting for both phonon and adsorbate (low-frequency 
modes) excitations at the same time is then easily obtained. Again, this theoretical 
procedure is more general than that previously reported based on phenomenological 
arguments [28]. 



H 




+ V(q) 




(34) 



N ph 2 N ad ^ 2 

2 niiuf 2 ^— ' 

i=i 1 7=1 j J 



(35) 
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(d) The convolution problem. Some simple applications 

As a final remark, we would like to mention that, taking into account all the 
previous information, one could carry out the calculation of the convolution of phonon 
and adsorbate lineshapes in order to extract reliable lifetimes and damping coefficients 
or phonon frictions in a more careful way. Thus, as an illustration, now we will 
focus on experimental lineshapes obtained from typical time-of-flight measurements 
from He atom-surface scattering: an isolated frustrated translational (T) mode peak 
and the masking of a T-mode by a substrate Rayleigh phonon. In particular, we will 
consider He scattering off a Cu(OOl) surface with a 3% of CO coverage at low (surface) 
temperatures [20], displaying the results in terms of the energy transferred (AE) in 
the scattering process. Experimentally, the intrinsic peak width is of the same order of 
magnitude as the apparatus energy resolution. Therefore, quantitative reliable results 
can only be extracted by using deconvolution techniques. Thus, in the case of low surface 
temperatures, one proceeds (numerically) by assigning a single Lorentzian function to 
the intrinsic lineshape and convoluting it with the apparatus response function, which 
is usually assumed to be a Gaussian function. Then, by means of standard least squares 
fitting techniques, the resulting curve is optimized until it fits the experimental one, 
extracting finally from this procedure the good Lorentzian lineshape. In the applications 
below, we have used this procedure rather than Eq. (133]) . since we are within a low 
temperature regime and only one Lorentzian function is going to contribute. Information 
about position of the maximum and full width at half maximum is processed in terms of 
renormalized frequencies according to (123|) and fl26|) . As the surface temperature, more 
and more Lorentzian functions have to be added in order to determine the lineshape, in 
agreement with (1331 . 

In Fig. [T] we show the result obtained (solid line) from the convolution of a single 
Lorentzian function with the Gaussian apparatus response function, i.e., a typical Voigt 
profile, which describes the creation of a T-mode peak at AE = hwo ~ —4 meV and a 
surface temperature T s = 60 K. The dashed line represents the best-fitting Lorentzian 
function whose position and full width at halh maximum (FWHM) are determined, as 
mentioned above, by convoluting the Lorentzian with the apparatus response function 
and then finding the best fit to the experimental data (here, they have been extracted 
from Ref. [20] and are represented with open circles). As can be seen, the agreement 
between the experimental data [20] and our fitting is fairly good. The FWHM is around 
0.38 meV (damping constant) and u can be calculated from (1231) . As mentioned above, 
as the surface temperature increases, the number of Lorentzian terms to be considered 
in (1331) will also increase in order to reproduce a wider and lower convolution profile, 
which can be later compared with the experimental peak. 

In Fig. El the results from the convolution (solid line) describe the case where 
the peaks corresponding to the annihilation of both a substrate Rayleigh phonon (at 
AE ~ 3.5 meV) and the T-mode of the adsorbate (at AE ~ 4 meV) are very close. 
Here, again, the surface temperature is T s = 60 K, which allows us to consider only one 
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Figure 1. Comparison between experimental data (o) [20] and convolution results 
(solid line) for the creation of the T-mode peak (AE w —4 meV) at T s = 60 K. 
The dashed line represents the best-fitting Lorentzian function whose position and 
FWHM are determined in such a way that, once the Lorentzian is convoluted with the 
(Gaussian) apparatus response function, it best fits the experimental data. 



Lorentzian function (dashed lines) for each event, the excitation of the substrate phonon 
and the adsorbate T-mode. As can be seen, we find a fairly good agreement between our 
theoretical results and the experimental data (open circles) [20]. Both peaks are so close 
that they overlap, which leads the Rayleigh phonon peak, more intense, to mask the 
peak corresponding to the T-mode. This is in a sharp contrast with the relative height 
of the corresponding Lorentzians, which is lower for the former. This overlapping of 
the two Lorentzian functions (dashed lines) thus explains the formation of the shoulder 
observed at larger AE. 
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